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Abstract 

Conventional weak-coupling Rayleigh-Schrodinger perturbation theory suf- 



^ ■ fers from problems that arise from resonant coupling of successive orders in 

■ the perturbation series. Multiple-scale analysis, a powerful and sophisticated 

f"*) . perturbative method that quantitatively analyzes characteristic physical be- 

0^ ■ haviors occurring on various length or time scales, avoids such problems by 

implicitly performing an infinite resummation of the conventional perturba- 
^ . tion series. Multiple-scale perturbation theory provides a good description 

of the classical anharmonic oscillator. Here, it is extended to study (1) the 
Heisenberg operator equations of motion and (2) the Schrodinger equation for 
the quantum anharmonic oscillator. In the former case, it leads to a system 
of coupled operator differential equations, which is solved exactly. The solu- 
tion provides an operator mass renormalization of the theory. In the latter 
case, multiple-scale analysis elucidates the connection between weak-coupling 
perturbative and semiclassical nonperturbative aspects of the wave function. 
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I. INTRODUCTION 



Multiple-scale perturbation theory (MSPT) is a powerful and sophisticated perturbative 
technique for solving physical problems having a small parameter e . This perturbation 
method is generally useful for both linear and nonlinear problems. In fact, it is so general that 
other well-known perturbative methods, such as WKB theory and boundary-layer theory, 
which are useful in more limited contexts, may be viewed as special cases of MSPT |]T| . 

The key idea underlying MSPT is that dynamical systems tend to exhibit distinct char- 
acteristic physical behaviors at different length or time scales. If a conventional perturbation 
series is used to solve a problem, then there is often a resonant coupling between successive 
orders of perturbation theory. This coupling gives rise to secular terms in the perturbation 
series (terms that grow rapidly as functions of the length or time variable). These secular 
terms conflict with physical requirements that the solution be finite. MSPT reorganizes 
the perturbation series to eliminate secular growth, and in doing so it gives a quantitative 
description of the characteristic behaviors that occur at many scales. In the past MSPT 
has been used to solve classical differential equations such as the equation of motion for the 
classical anharmonic oscillator. Indeed, the classical anharmonic oscillator is often used to 
illustrate and explain the method of MSPT. 

In this paper we generalize the ideas of MSPT to the quantum anharmonic oscillator . 
The quantum anharmonic oscillator is an excellent laboratory for the study of a variety of 
perturbative methods. It has been used to study the origin of the divergence of conventional 
weak-coupling Rayleigh-Schrodinger perturbation theory Pade and Borel summation of 
perturbation series 0, large-order behavior of perturbation theory [0, delta expansions [§, 
dimensional expansions P], and strong-coupling expansions [jlO|- Here, we use MSPT to 
study two aspects of the quantum anharmonic oscillator, the Heisenberg operator equations 
of motion in Sees. |I| and |T| and the Schrodinger equation in Sees. |^ and 0. We illus- 
trate the methods of MSPT in Sec. |I| by applying it to the nonlinear dynamical equation 
of motion for the classical anharmonic oscillator (Buffing's equation). There, we obtain 



2 



the first-order frequency sliift. Tlien, in Sec. |TT| we extend tlie metliods of MSPT to the 
nonlinear Heisenberg equation for the quantum anharmonic oscillator (the quantum version 
of Dufhng's equation). To complete the analysis it is necessary to solve a nonlinear system 
of coupled operator differential equations. We find the exact closed-form solution to this 
system. From this solution, we obtain the quantum operator analog of the frequency shift; 
namely, an operator mass renormalization that expresses the first-order shift of all energy 
levels. 

In the next two sections we study the Schrodinger equation for the quantum anharmonic 
oscillator. Specifically, we examine the asymptotic behavior of the wave function tplx) for 
large x. We consider the problem of reconciling the different results that one obtains from 
conventional Rayleigh-Schrodinger perturbation theory (a formal Taylor series in powers of a 
small parameter e) and WKB theory (a nonperturbative probe of the anharmonic oscillator 
that is valid regardless of the size of e). To any finite order in conventional perturbation 
theory, tplx) behaves like a Gaussian for large x; however, WKB theory predicts that as 
X — oo the wave function decays to zero like the exponential of a cubic. In Sec. |V| we 
resolve this discrepancy at two different length scales by an infinite sequence of reorderings 
and resummations of the conventional weak-coupling perturbation series. In Sec. [V| we 
explain the origin of the disparity by performing a direct multiple-scale analysis of the 
Schrodinger equation for the quantum anharmonic oscillator. 

The approach used in this paper for the anharmonic oscillator wave function has been 
applied in perturbative quantum field theory to sum leading-logarithm divergences and 
leading infrared divergences |jl2|]. It is our hope that in the future the direct nonperturbative 
multivariate approach of MSPT will provide a framework to simplify such schemes. 
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II. MULTIPLE-SCALE PERTURBATION THEORY APPLIED TO THE 
CLASSICAL ANHARMONIC OSCILLATOR 



In this section we explain MSPT by using it to treat the classical anharmonic oscillator, 
a dynamical system satisfying the nonlinear differential equation 



which is known as Buffing's equation. The positivity of e ensures that there are no un- 
bounded runaway modes. We impose the initial conditions 



The harmonic oscillator (e = 0) has only one time scale, namely, the period of oscillation. 
However, the nonlinear term in Eq. (|2.1|) introduces many time scales into the problem. For 
example, when e 7^ 0, one can observe on a long-time scale [t = 0(e"^)] a frequency shift 
of order e. One can study the classical anharmonic oscillator on the short-time scale t and 
also on many long-time scales r = et, Ti = e^t, T2 = e^t, and so on. 

Let us first examine what happens if we attempt to solve Buffing's equation using a 
conventional perturbation series in powers of the parameter e. 




(2.1) 



y{0) = 1 and y'{0) = 0. 



(2.2) 



00 



(2.3) 



n=0 



in which the initial conditions in Eq. ( p.2|) are contained as 



yo(0) 



1 and y'o{0) = 0, 




(2.4) 



Substitute Eq. (pj) into Eq. (gjj) 



To leading order (zeroth order in powers of e), we have 




(2.5) 



and to first order in powers of e we have 
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^1/1 + y, = -Ayl (2.6) 



The solution to Eq. ( |2.5| ) satisfying the initial conditions in Eq. ( p.4| ) is 

2/o(t) = cost (2.7) 

When we introduce this solution into Eq. (|2.6|) , we obtain 

(f 

^Z/i + = -cos(3t) -3cost. (2.8) 

Equation ( |2.8D represents a forced harmonic oscillator whose driving term has frequencies 3 
and 1. When a harmonic oscillator is driven at its natural frequency, which in this case is 
1, we have the phenomenon of resonance. As a result, the solution 

1 13 

?/i(t) = - cos(3t) cost tsint (2.9) 

8 8 2 



contains a secular term that grows linearly with increasing time t. Equation ( |2.9| ) cannot 
be valid for long times because the exact solution to Buffing's equation remains bounded 
for alH |]1[]. Hence, the conventional perturbation expansion is sensible only for short times 
t « e^^. How then does the conventional perturbation series determine the behavior of 
y{t) for long times, say of order e~^? 

One way to answer this question is to identify the structure of the most secular (highest 
power in t) term to all orders in perturbation theory. One can easily verify [|l| that for all n 
the most secular term in ?/„(t) has the form 



ir 

2^ 



3A" , / 3^\" 



27 ^ n-T' ' 



(2.10) 

Since the expression in Eq. (|2.1CI|) is multiplied by e", if we make the approximation of 



retaining only the most secular term in every order, then we obtain a series in powers of the 
long-time variable r = et. Evidently, we can sum the most secular terms in this series to all 
orders in e, and since the result is a cosine function, it remains bounded for all times t: 



We interpret this result to mean that on the long-time scale r there is a frequency shift in 
the oscillator of order |e. Of course, this result is not exact; there are less secular terms 
to all orders in the perturbation expansion, and these terms give rise to frequency shifts of 
order e^, e^, and so on. 

We will now show how MSPT directly reproduces the result in Eq. ( |2.11| ). To avoid 



the complicated procedure of summing the conventional perturbation series to all orders 
in powers of e, MSPT uses a sophisticated perturbative approach that prevents secular 
terms from appearing in the perturbation expansion. Multiple-scale analysis assumes a 
priori the existence of many time scales {t, r, ri, T2, ...) in the problem, which can be 
temporarily treated as independent variables. Here, we illustrate by performing just a first- 
order calculation. We use only the two variables t and t = et and seek a perturbative 
solution to Eq. ( |2.1j ) of the form 

y{t) = Yo{t, t) + eY,{t, r) + 0{e^). (2.12) 

Using the chain rule and the identity % = e, we convert Eq. to a sequence of partial 



at 

differential equations for the dependent variables Yq, Yi, ... . The first two equations read 

^n + io^O, (2.13) 
^y-+yi._«--2^y, (2.14) 



The general solution to Eq. ( p.l3| ) is 



Yo{t,T) = A(r) cost + 5(r) sin t, (2.15) 

where A{t) and B{t) are as yet unknown functions of r. We determine these functions by 
imposing the condition that no secular terms appear in Yi. That is, we use the functional 
freedom in the choice of A{t) and B{t) to eliminate the resonant coupling between zeroth 
and first order in perturbation theory. We substitute YQ{t,T) in Eq. ( |2.15D into the right 
side of Eq. ( p. 141) and expand the resulting expression using the trigonometric identities 
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1 3 

cos^ t = - cos(3t) + - cos t, 

4: 4: 

sm t = —- smlSt) + - smt, 

cos t sin t = - sin( 3t) + - sin t, 
4 ^ ^ 4 

sin^t cost = — - cos(3t) + - cost. (2.16) 

To eliminate secularity we require that the coefficient of cost and sint vanish. This gives 
the pair of equations 

2^ = -3A^-3AB^ and 2^ = 3B^ + SA^B. (2.17) 
dr dr 

To solve this system we multiply the first equation by B{t) and the second equation by 
A{t). Adding the resulting equations gives 

^C(r) = 0, (2.18) 

where 

C{r) = l[A{r)f + \[B{r)]\ (2.19) 

Hence C(r) is independent of r and we may take C(r) = C(0). 

Substituting this result back into Eq. ( |2.17D gives the elementary linear system 



^B = -3C{0)A and ^A = 3C{0)B. (2.20) 
dr dr 

When we solve this system and then impose the initial conditions, we obtain C(0) = | and 



Fo(t, t) = cos 



(2.21) 



where we have used r = et. We have thus reproduced the approximate solution in Eq.( p.llD 



that is valid for long times. Note that while conventional perturbation theory is only valid for 
t << e^^, the multiple-scale solution is valid for times satisfying t << e^^. This multivariate 
approach can, in principle, be performed to nth order in powers of e for any n. 
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III. MULTIPLE-SCALE PERTURBATION THEORY APPLIED TO THE 
QUANTUM ANHARMONIC OSCILLATOR 



In this section we extend the multiple-scale approach described in Sec. y to the quan- 
tum anharmonic oscillator. This is a nontrivial generalization of the usual multiple-scale 
techniques because it requires that we solve operator differential equations ||1^ . 



The quantum anharmonic oscillator is defined by the Hamiltonian \T^\ 



H{v.q) = \f + \q' + ^q' (e > 0), (3.1) 

where p and q are operators satisfying the canonical equal-time commutation relation 

[q{t),p{t)]=in. (3.2) 

The positivity of e ensures that H{p, q) is bounded below. 
The Heisenberg operator equations of motion are 

|,= i-[g,i7(p,g)]=p, (3.3) 
j^p = ^[p,H{p,q)] = -g(t) -4e[g(t)]l (3.4) 

These equations combine to give 

^q{t)+qit)+Ae[qit)]' = 0, (3.5) 



which is the quantum analog of Buffing's classical differential equation ( p.l|) . Since p{t) and 
q{t) are operators, we cannot impose numerical initial conditions like those in Eq. ( p. 21) ; 
rather, we enforce a general operator initial condition at t = 0: 

g(0) = go and p{0) = po- (3.6) 

Here, po and go ai"e fundamental time-independent operators obeying the Heisenberg algebra 

[qo,Po\ = ih. (3.7) 
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Rather than seeking a solution to Eq. (p.5|) as a conventional perturbation series in 
powers of e, we perform a multiple-scale analysis. We assume that q{t) exhibits characteristic 
behavior on the short-time scale t and on the long-time scale r = et and we write 

g(t) = g(t, r) = Qo(t, r) + eQ^it, r) + Oie^). (3.8) 

This equation is analogous to Eq. (|2.12| ) but here Qo and Qi are operator-valued functions. 
There is an associated expression for the momentum operator p{t): 

p{t) = Po{t,T) + eP,{t,T) + 0(e2). (3.9) 

Furthermore, since the momentum operator is the time derivative of the position operator 
[see Eq. ( p.3|) ], the chain rule gives 

p{t) = ^t^o + e (-^Qo + ^Qi) + 0{e'). (3.10) 



We substitute the expression for q{t) in Eq. ( |3.8| ) into Eq. (|3.5|), collect the coefficients 
of e*^ and e^, and obtain operator differential equations analogous to Eqs. ( |2.13| ) and ( ^.14^ : 



Q2 

_go + go = 0, (3.11) 

§^Qi + Qi = -^Ql-2^Qo. (3.12) 

Because Eq. (|3.11|) is linear it is easy to find its general solution: 

Qo{t,T) = A{t) cost + B{t) sint; (3.13) 

from this result we obtain 

Po(t,r) = i3(r) cost -^(r) sint. (3.14) 

It is now necessary to find the coefficient functions A{t) and B{t), which are operators. 
The canonical commutation relations in Eq. (|3.2[) imply that these operators must satisfy 



[Air),B{T)]=ih. (3.15) 



Also, the initial conditions in Eq. 



give 



^(0) = go and B{0) = po. 



(3.16) 



To determine the r dependence of the functions A{t) and B{t), we must eliminate secular 
behavior. To do so we examine the right side of Eq. ( |3.12| ), which is 



4 [A{t) cos t + B{t) sin - 2 



--^A{t) sin it + -^B{t) cost 
ar dr 



(3.17) 



Next, we expand the cubic term, taking care to preserve the order of operator multiplication, 
and we use the trigonometric identities in Eq. (|2.16|) . We eliminate secularity by setting the 
coefficients of cost and sint to zero and obtain 

dB dA 
2— = ~3A^ - BAB - BBA - ABB and 2— = 3B^ + ABA + AAB + BAA. (3.18) 
dr dr 

This system of operator- valued differential equations is the quantum analog of Eq. ( |2.17| ). 



To solve the system ( p.l^ ) we begin by multiplying the first equation on the left and on 
the right by B{t) and the second equation on the left and on the right by A{t). Adding the 
resulting four equations and simplifying, we get 

d 



dr 



n = o, 



(3.19) 



where 



n^\A' + \b\ 



(3.20) 



Equation (p.l9|) is the quantum analog of ^ = in Eq. ( |2.18|) . By construction, the 
operator is independent of the short-time variable t. However, Eq. ( p.l9| ) shows that li. 
is also independent of the long-time variable r. Therefore, Eq. ( |3.16| ) allows us to express 
a. in terms of the fundamental operators po and q^: 



^ = 2^0 + 2^0- 



(3.21) 



Next we use Eq. ( |3.15|) to rewrite Eq. (|3.18|) in manifestly Hermitian form: 
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7 Q Q 

B = --{HA + AH) and —A=-{nB + Bn). (3.22) 
CLT 2 dr 2 

Suppose for a moment that we could replace Ti by the numerical constant C(0). We would 

then obtain the elementary classical system of coupled differential equations in Eq. ( p. 20] ). 

Because this system is linear we could treat these operator differential equations classically 

and ignore operator ordering. The solution to this system that satisfies the initial conditions 

in Eq. ( |3.16| ) would then be 

A{t) = go cos[3C(0)r] + po sin[3C(0)r] and B{t) = po cos[3C(0)r] - go sin[3C(0)r]. (3.23) 

This solution suggests the structure of the exact solution to the operator differential 
equation system (|3.22|) . The formal solution is a natural quantum operator generalization 
of Eq. (|3.23|) using Weyl-ordered products of operators: 

A{r) = W[qo cos{3nT)] + W[po sin(37^r)], 

B{t) =W[po cos(3Hr)] - >V[go sin(3Hr)] . (3.24) 

The notation W[go/(7ir)] represents an operator ordering defined as follows: First, expand 
the function fiHr) as a Taylor series in powers of the operator Ht. Then Weyl-order the 
Taylor series term-by-term: 

W{qonn = ^ [(o)go^" + (i)^go^""' + (2)^'go^"-' + ... + C)^"go] . (3.25) 

Using this definition it is straightforward to verify that Eq. ( |3.24| ) is indeed the exact operator 
solution to Eq. ( |3.22| ) satisfying the initial conditions in Eq. ( p.l6| ). 



Our objective is now to simplify the formal solution in Eq. ( p.24| ) and to re-express it 
in closed form. Since sines and cosines are linear combinations of exponential functions, we 
consider first the general problem of simplifying the Weyl-ordered product 

1 ... 1 



Wfgoe^n = W 



(3.26) 



For each power of r we reorder the operators by commuting go symmetrically to the left and 
to the right to maintain the Hermitian form 
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>V(go) = go = 2(90 + Qo), 



W(goH) 



2 v^°¥ 



■go 











(re 




' 2 






-a- 




- i) 



2 



go 



IT 



"W 



IT 



2 



2 ^2 



16 j + I ^4 



2^2 



16 



go 



and so on. This process defines a set of polynomials in the variable Ti/h ||I5| 



1, 






u 








1 






~ 4' 
3H 




' 


3 7^2 5 

~ 27F^ 16' 
5 25 n 




Y' 




re 


~ 2^^ 16¥' 

15 75 


61 




64 


Is as 


Euler polynomials 


m 



(3.27) 



(3.28) 



^: En + l)- The generating function for these nonorthogonal polynomials is given by 

^-i:^^»(^i) (M<-). (3.29) 
This generating function allows us to express the following Weyl-ordered product compactly: 



W goe 



(3.30) 



2cosh(r;i/2) 

Using Eq. ( p. 30] ) the cosines and sines of our quantum solution in Eq. ( p. 24] ) can also be 
written in compact form when we take combinations of complex exponentials: 
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W gocos 37^r = , , , , 

2 cos(3rn,/2) 

^A)l ■ Ci'u M gosin(37^r) + sin(37^r)go 

W gosin 3Hr = , , , . 3.31 

2 cos{6Th/2) 

Last, we substitute this result into the zeroth-order solution in Eq. ( ^.13| ) and obtain 

, , _ go cos{t + SHt) + cos(t + 37Yr)go Po sin(t + 3Ht) + sm{t + 3'Ht)pq 
' " 2cos(3r;i/2) ^ 2cos(3r;i/2) 

go cos(t + 3Het) + cos(t + 3Het)qQ pq sm(t + 3Het) + sin(t + 3Het)pQ 
~ 2cos(3et^/2) ^ 2 cos(3et;i/2) ' ^' ^ 

where we have replaced r by et. The result in Eq. ( |3.32|) is the objective of our multiple-scale 
analysis. It is the quantum operator analog of Eq. (|2.21| ). Indeed, in the limit as h ^ 0, we 
recover the classical multiple-scale approximation in Eq. ( p.21|) . [To obtain this result we 
impose the classical initial conditions po = and go = 1, which from Eq. (|3.2|) give Ti. = |.] 

We interpret this multiple- scale approximation to the operator g as follows: In the clas- 
sical case we identify the coefficient of the time variable t as a first-order approximation to 
the frequency shift. However, here the coefficient of t in Eq. ( |3.32| ) is an operator. Thus, we 
have derived an operator form of mass renormalization. 

To understand this operator mass renormalization we must take the expectation value of 
Eq. ( |3.32| ) between states. Then by examining the time dependence of this matrix element 
we can read off the energy-level differences of the quantum system. It is easy to construct a 
set of states because the operators go and po satisfy the commutation relation (|3.7|) . Hence, 
appropriate linear combinations of go and po may be used as raising and lowering operators 
to generate a Fock space consisting of the states \n). By construction, these states are 
eigenstates of the operator H: 

n\n) = (n + ^^h\n). (3.33) 

Let us take the expectation value of Eq. ( |3.32| ) between the states {n — 1\ and \n). 
Allowing the operator Ti, to act to the left and the right using Eq. ( p. 33 ), we obtain 

t + 3(n+fjhet +COS t + 3 (n - |) het 



cos 

{n - l\Qo\n) = {n- l|go|n) 



2cos(3et/i/2) 
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+ {n - l\pQ\n)- 



sm 



t + ?,(n + fjhet + sin t + 3 - ^et 



2cos(3et;i/2) 

{n — llgol'^) cos[t(l + ?>nhe)] + (n — l|po|^) sin[t(l + 3n^e)], 



(3.34) 



and we can see that the energy-level differences of the quantum oscillator are 1 + "inhe. 

Let us verify this result. If we set e = in Eq. p.l| ), we obtain the harmonic oscillator, 
whose coordinate-space eigenfunctions ipn{x) and corresponding energy eigenvalues En are 



_1 2 1 

i^nix) = e 4^ He„(x) and -E^ = n + -, 



(3.35) 



where He„(x) represents the Hermite polynomial. The first-order perturbative correction to 
the energy eigenvalues is obtained by computing the expectation value of the perturbation 
term in the Hamiltonian H: 

1 e;i/°°,rfxe-5^'[He„(x)]V 1 3,, , x 2x . x 

En = n + - + .LJ ... =n + - + -eh{2n^ + 2n + 1) + 0{e^). (3.36) 



2 4 /-^rfxe-2-^[He„(x)] 



2 4 



If we now calculate the energy difference En — En-i from Eq. (|3.36| ), we obtain 



En - En-i = 1 + 3nhe + 0(e2), 



(3.37) 



which verifies the result in Eq. (p.34|) . 



IV. RESUMMATION OF THE CONVENTIONAL PERTURBATION SERIES 
FOR THE QUANTUM ANHARMONIC OSCILLATOR 



Now we turn our attention from the Heisenberg equations of motion to the Schrodinger 
equation and examine the behavior of the wave function. The ground-state wave function 
tplx) for the quantum anharmonic oscillator satisfies the Schrodinger equation 



(4.1) 



and obeys the boundary conditions 



'4){±oo) = 0. 
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(4.2) 



We can use WKB analysis to find the large-x asymptotic behavior of the wave func- 
tion iIj{x). a geometrical-optics approximation gives the controlling factor (the exponential 
component of the leading asymptotic behavior) as 

e-v^l-lV6_ (4,3) 

This result is nonperturbative in the sense that WKB is valid independent of the size of the 
parameter e. The question addressed in this section is whether it is possible to reproduce 
this result using perturbation theory. 

The conventional approach to solving Eq. ( [4.1|) using Rayleigh-Schrodinger perturbation 
theory |^,|^ represents both the eigenfunction and eigenvalue as power series in e: 

oo oo 

^(x) ~ ^ e>„(a;) and E(e) ~ J] e"E„. (4.4) 

n=0 n=0 

Note that we use the asymptotic symbol ~ because, as is well known, the conventional 
perturbation series for the anharmonic oscillator diverges. 

Equations ( ^.Ij ) and ( ^.2| ) are homogeneous, so we are free to adopt the normalization 
condition ip{0) = 1, which translates into 

l/o(0) = l and |/„(0) = {n > 0). (4.5) 

The unperturbed (harmonic-oscillator) solutions corresponding to e = are 

l/o(x) = e-^'/^ and Eo = ^. (4.6) 

In Ref. it is shown that for all n, yn{x) is a product of the zeroth-order approximation 
to the wave function Uq^x) given in Eq. (^^ ) multiplied by a polynomial Pn{x): 

= e-^'/^P„(x). (4.7) 

Substituting Eq. (f4.7|) into Eq. (^?T|), we obtain the recursion formula for the polynomials: 

-1 n— 1 

P:{x) - xP'^x) = -x'Pn-i{x) - Y: PA^)En-r (4.8) 
^ i=o 
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The form of this recursion relation is generic. It is a typical recursive structure that arises 
in all perturbative calculations; to wit, the homogeneous part of this recursion relation is 
the same for all n while the inhomogeneous part contains all previous polynomials. It is this 
kind of recursive structure that is responsible for successive orders of perturbation theory 
being resonantly coupled. Here, the resonant coupling causes the degree of the polynomials 
to grow with n\ Pn{x) is a polynomial of degree 2n in the variable x"^: 



Po{x) 
Pi{x) 

P2{X) 

Psix) 
Piix) 

P,(x) 



3 /xN 2 / X 

2 [2) " [2 

31 13 1 /X 



21 /x^ 2 



+ 



6 



+ 



2 V2 



333 243 271 /x^ ^ 

,2 



30885 /xV 



+ 



64 
7 



2I + 



2777 ^ 



12 \2 



X 



14 



+ 



24 



16 

18461 



47 

2, 



17 fx^ 10 
12 V2 



X 



96 



2I + 



9195 /X 
128 V2 



1 

6 V2 

979 fx-- 10 



12 



48 



599 fx 



916731 651363 fx\^ 



128 
29177 



384 



128 
1325 /x^ 14 

~96~ U 



89673 

~ 32 V2/ 
269 25 



69107 fx^^ 
,2. 



64 



144 V2 



250183 fx^ 10 

,2. 



12 



144 V2 



144 V2 



X 



18 



1 

120 



768 

20 



Pdx) 



65518401 23046319 fx\^ 75770813 fx^ ^ 2476011 fx^^^ 



+ 



512 V2; ^ 
9259481 fx^ 10 



1536 
4055 /X 



256 V2; ^ 
13796435 fx^ 12 



+ 



648 



18 



+ 



^ 9216 
349 fx^ 20 



576 V2 



2 

29 



+ 



1536 
77173 fxV^ 



X 

720 V2 



256 

22 ]^ 



X 



For n > the general form of the polynomial is 

2n 



Pn{x) = Y.Cn,k{-: 



■X 



720 V2 



2\k 



+ 



24 



128 

112483 fx^ 16 



2304 



(4.9) 



(4.10) 



k=l 



Note that by virtue of Eq. (|4.5| ) the polynomials Pn{x) have no constant term when n > 0. 

We can derive a formula for En and a recursion relation for the coefficients Cn,k by 
substituting Eq. ([4.10|) into Eq. (|48D to obtain for n > 

Cn,l = Enj 

n-1 

2kCn,k + Cn-i,k-2 = -{k + 1)(2A; + l)Cn,k+i + E C„fca-,,i. (4.11) 
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The coefficients Cn,k form a triangular array in the sense that the degree of the polyno- 
mials increases with n. The convolution in Eq. ( |4.11| ) makes this recursion relation highly 



nontrivial. However, for all n it is possible to find C„.2n, the coefficient of the highest power 
of X. This is an exact analog of finding the coefficient of the highest power of t (most secular 
term) in nth order in perturbation theory for the classical anharmonic oscillator. By setting 
k = 2n in Eq. (|4.11| ), we see that Cn,2n satisfies the simple linear recursion relation 

4nC„,2„ + Cn-l,2n-2 = 0, (4. 12) 

whose solution is 

/ 1\" 1 

C,.a„ = (-j) (4,13) 

To study the behavior of the wave function ip{x) for large x, we approximate ipi^) by 
resumming the perturbation series in Eq. ([4.4|) and keeping just the highest power of x in 
every order. Using Eq. ( [4.13| ) we obtain a simple exponential approximation to tp{x): 



n=o n.lG" 

For the classical anharmonic oscillator this approach also gives an exponential approximation 
[see Eq. (|2.11|) ]. However, the classical and quantum anharmonic oscillators are evidently 



quite different; although we have summed the most secular terms to all orders in perturbation 
theory, the result in Eq. ( [4.14|) is not the correct behavior of the wave function ip{x) for 



large x because it decays to zero too rapidly. The correct behavior from WKB theory is an 
exponential of a cubic [see Eq. (|4.3| )] and not an exponential of a quartic! 

Can we improve our estimate of ifj^x) by including C„,2n-i, the coefficient of the next 
highest power oi x, in our summation? This coefficient satisfies the inhomogeneous difference 
equation obtained by setting k = 2n — 1 in Eq. ( [4.11| ): 

/ 1\" 1 

(4n - 2)C„,2„„i + C„_i,2n-3 = -2n{An - 1)C„,2„ = -2n{An - 1) [--j -. (4.15) 

The solution to this inhomogeneous equation is 
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IX"™ 1 n 

Cn,2n-l = -[—^) ^^3(4^ + 5). (4.16) 



If we include this formula in the resummation to all orders, we obtain the result in Eq. ( [4.14| ) 
now multiplied by a polynomial: 



1 



We have again failed to obtain the correct large-x behavior of ip{x). Nevertheless, let us 
continue to reorganize the perturbation series. Setting k = 2n — 2 m Eq. (f4.11|) , we get 



3 

{An - A)Cn,2n-2 + C„_i,2n-4 = -(2n - l)(4n - 3)C„ 2„-1 + -C„_i,2n-2, (4.18) 



1^" 1 n(n- 1),^^ 2 



whose solution is 

C„,2n-2 = (^--j ^ '"^'"^g ^' {lQn' + Un + %l). (4.19) 

The next few coefficients are: 

/ l\"lr?fr? — l") 
C„,2n-3 = - -7 -7 ^ (64n^ + 400^3 + 764^^ - 433n + 390), 

V 4/ n! 162 



n,2n-4 



1\" 1 n{n-l){n-2) 



n,27i— 5 



4/ n! 1944 
X (256/2^ + 2816n^ + 11744n^ + 18304^2 + 34209n + 70029), 
1\" 1 n{n - l){n - 2) 



4/ n\ 29160 
x(1024n^ + 14592n^ + 74368n^ + 106080n^ + 45316n^ 

+143073^2 - 2392803n + 3967380), 
1\" 1 n(n- l)(n-2)(n-3) 



<^-.2n-6 ^ 4; ^, 524880 

x(4096n^ + 86016n^ + 718080n^ + 2799360ra^ + 6702384n^ 

+16486704^3 + 16745975n2 + 180087585n + 415966860). (4.20) 

Evidently, the coefficient Cn,2n-j has the general form of a polynomial in n of degree 2j 
multiplied by (— 4)~"/?t,!. 

When we sum these coefficients to all orders, we obtain a totally new sequence of poly- 
nomials (displayed below in square brackets) in the variable ex^: 



2 V2/ 3 



+ 
+ 
+ 

+ 

+ 



■31 9 /xX 

r (2) 

21 , /x\ 

,r (2) 

243^3 



2 2 /x^ 



2 4 /X' 



,4 5307 4 
' ^ 128"^ 1 - ) - 

, 20 



12 



10 



14 



18 



27 V2/ 243 V2. 
333^3 (xy ^ 14465^4 /x' 



24 



133 6 
'36^ V2y 



16 



39493 5 /xyo 12463^ 
256 ^ V2/ ^ 192 



313 7 fx^ 



18 



96 



+ 



211 8 /X' 



162 



■2777 4 /xy 46891 5 
" ^ 2/ 6^^ V2 



X 



.2, 

, 8 



X 



14 



5444579^6 



64 " K2J ^ 9216 
x\20 833 q /x\24 307 ^ /2;\28 

'2) + 972' (2) 



x\i2 14497 7 
e 



79357 8 /a;\20 833 

^ 1728 ^ V2/ ~ 162' 

>32 4 (^X 



64 



30 



X 



16 



4 ^1 /X\32 4 

"405^ U/ ^ 32805^ 



, 36 



(4.21) 



The original form of the perturbation series in Eqs. ( |4.4| ) and ( [4.7[ ) has undergone a re- 
markable transmutation. The original polynomials Pn{x) have been absorbed and the wave 
function iIj{x) is now represented as a more elaborate exponential multiplying a new class of 
polynomials. [We can, of course, generate these polynomials directly from the Schrodinger 
equation (|4.1| ) using a recursion relation similar to that in Eq. (f4.8|) .] 

Let us perform a further resummation procedure. That is, for this new set of polynomials 
we determine the coefficient of the highest power of x, 

'2\" /x\6" 1 



37 V2y n!^ 



the coefficient of the second-highest power of x, 

9 /2\" /x\6"-4 1 



2n 



(4.22) 



4 V3 



n! 



(4.23) 



the coefficient of the third-highest power of x, 

I) 



9 /2 
32 V3 



ry, \ 6n— 8 1 

^ -e'"-2n(n-l)(9n2-3n+l), 
n! 



(4.24) 



19 



the coefficient of tfie fourtfi-fiigfiest power of x, 

"128 (3) (2) ^e'--Mn-l){n-2){9n'-9n' + 7n + 4), (4.25) 
and so on. 

Having discovered tfiese formulas we can now sum over n. This further reorganization 

of the perturbation series gives a new approximation to iIj{x) as an exponential multiplied 
by yet another set of polynomials, this time in the variable e^x^: 



g-x2/4g-ex4/16ge2a;6/96^ ^ _ 

+ 



3 / n f X 



31 2 , 7 „ 



X 

2 V2 



10 



-e , 
2 V2 



16 



187 o fx\^ 277 . fx 



12 



11 7 fxy^ 1 „ 

(2) +r (2 



24 n 



+ ... . 



(4.26) 



For these new polynomials the term containing the highest power of x has the form 



(4.27) 



and the term containing the second-highest power of x has the form 

8n-6 



-(-irf- 

2^ ' V2 



e^"-^n(4n- 1). 



(4.28) 



Again, we reorganize the perturbation series by summing over all n. This resummation 

gives yet another representation for il){x) as an exponential multiplied by a new set of 
polynomials, this time in the variable e^x^: 



g -a;2/4g-ea;'*/16ge2a;S/96g-€3a;8/256 J 



10 



+ 



2 \2J \2, 

For these polynomials we identify the term containing the highest power of x. 



(4.29) 



lOn 



2" - — e*". 



(4.30) 



and we reorganize the perturbation series again by summing over all n. This resummation 
gives a new exponential and a new set of polynomials, this time in the variable e^x^^: 
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g-a:2/4g-ex'l/16ge2a;6/96^_e32:8 



/256ge4.iO/5i2|^^ _ j_ (4.31) 



We can continue this process indefinitely. With each iterative reorganization of the 
perturbation series we generate a new exponential multiplied by new polynomials. However, 
we have not attained our original goal of deriving the exponential behavior in Eq. ( [1.3[ ) from 
the weak-coupling perturbation series in powers of e. Indeed, it seems impossible for our 
approach to succeed because at each stage in the reorganization of the perturbation series 
the variable x appears only in even powers. How can we ever obtain the exponential of a 
cubic! 

There is a simple and direct answer to this question. We merely recognize that the 
exponent in Eq. ( [4.31| ) is the beginning of a binomial series: 

X 1 /( Itr l^is l/iin 1 



1 - 1 + ex 



2\3/2 



4 16 96 256 512 6e 

If we now let x be large (ex^ >> 1), we recover the asymptotic behavior in Eq. 

1 



(4.32) 



exp 



6e 



1 - 1 + ex 



2\3/2 



Let us examine more deeply what happens at each reorganization of the perturbation 
series. At the first resummation we sum over terms of the form e"x'^" and treat ex^ as small. 
Thus, while x may be large compared with 1, it cannot be too large; x must satisfy the 
asymptotic bound x << e~^/^. At the next reorganization we sum over terms of the form 
^2n^6n treat e^x^ as small. This resummation is valid in the larger region x << e^^/^. 
At the next level we sum over terms of the form e^^x^'"' and treat e^x^ as small. Hence, this 
resummation is valid in a still larger region, x << e~^/^. At the jth iteration the range of x 

j 

increases to x << e^j+a . Clearly, as j oo, we obtain an estimate of the wave function ipi^x) 
that is valid for x as large as e^^^^. It is only when x is this large that the ex^ term in the 
Schrodinger equation ( [4.1| ) becomes comparable in size to the x^ term! Thus, after a finite 



number of reorderings of the perturbation series, we cannot expect to reproduce exactly the 
exponential behavior in Eq. 
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Nevertheless, at each stage of the resummation process we observe precursors of the 
exponential behavior in Eq. ( [4.3| ). To demonstrate this, we examine the structure of the 
exponent in Eq. ( [4.31[ ), in which we factor out the first term: 



+ iex^ - -e'x' + -eV - — + ..) . (4.34) 
4 V 4 24 64 128 / ^ ^ 

At the first iteration, which is valid for x << e"^/^, we neglect all terms in Eq. ([4.34|) beyond 
^ex^ (because they are small compared with 1) and replace this series by 

^2 . 1 .1/2 

-^(l + M . (4.35) 



At the second iteration, which is valid for x << e we neglect all terms in Eq. (|4.34|) 



beyond —^e'^x^ and replace the series by 



4 V 24 
At the next two iterations we have 



^ + ea;2^AeV) . (4.36) 



X 



11 Q \ 1/6 

-L-LO/1 ^ Q a \ 



^ ( 1 + ^ex' + ^e'x' + ^eV ) (4.37) 



and 



-^(l + 2e.= + ^.V + A,3,.^JZ_,V)"'. (4.38) 

Now, let X be large. In each of the above formulas we obtain a cubic in x\ Moreover, 
this cubic is multiplied by —y/e and the numerical coefficients of the cubic approach |: 

1 1 



4^2 5.65685' 
1 1 



4(24/5)1/4 5.92066' 
1 1 



4(32/3)1/6 5.93469' 

^ ^ (4.39) 



4(1152/47)1/8 5.96663' 
It is now clear how it is possible to obtain a cubic in x from only even powers of x 
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Having given a detailed discussion of the exponential prefactor that emerges at each 
stage of the resummation process, we conclude this section by describing the structures of 
the polynomials that arise at every stage. Let N represent the number of resummations 
that we have performed. Then for each value of N we can express the wave function ip{x) 
as an exponential multiplied by a sum over polynomials. We have established that after N 
resummations the exponential factor has the form 



exp 



^2 N 



Hi) 



2\n 



r 1 3 

n=0 -L 



{ex') 



(4.40) 



2 nj(n + l)! 

At each stage N in the resummation process the argument of the polynomials changes so 
we denote the argument by z^, where 

2\ A'+l 



(4.41) 



In general, the wave function is a product of the exponential factor multiplied by a power 
series in the variable 



_1 2 

y = -ex 



(4.42) 



We represent the coefficient of at the A^th resummation by the notation V^^\zn). Thus, 



ip{x) — exp 



•2 N 



[ex 



„2\n 



(4.43) 



n=0 



4 ^or(l-n) (n+1)! 

We now describe the structure of 'P^\zn)- For sufficiently large N, N > n — 2, 
is a polynomial of degree n. However, when N < n — 2, negative powers are present. The 
largest negative power of zn is zjl''^~^~^\ Thus, in the range from A^ = OtoA^ = n — 2, 
V^\zn) contains all powers of zn from 

^-{n-N-l) ^„ _ ^^^^ N>n-2, Vi^\zN) 

contains all powers of 2;jv from z% to z^. 

We now give some formulas for the functions P^^^ . For all A", v'q'^ is extremely simple: 



(4.44) 



We have found a closed-form general expression for V^'' valid for all A^: 
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-2{-l 



{2N)\ 



N\ {N + 2)\ 



(4.45) 



The closed-form expression for Pf ) is more complicated: 



2^0 



31 

8 



.NjV (lliV+13)(2jV)! ^ 



{2Ny. 



N\(N+2y. 



(iV = 0), 

4 (A^>o). 



(4.46) 



The expression for Vf' is still more complicated: 





333 ^-2 243 ^-1 271 
8 ^0 8 ^0 16 


47^ 


17^2 
12^0 


1 ^3 
6^0 


{N = 0), 






21 -1 187 _|_ 73 ^ 
4 ^1 16 12^1 


1 ^ 81 1 




{N = l), 


(4.47) 




~W + ^^^^ + '^^^Af ~ (~ 


1)-| 


(2Af)! 




(iV>l), 






Ar!(Ar+2)! 





where and cat are rational numbers. 



Ginsburg and Montroll [0 attempted to approximate the wave function for the anhar- 
monic oscillator by a sequence of exponentials of fractional powers of polynomials. In this 
respect their work is superficially similar to ours. However, they did not generate these 
polynomials from the Rayleigh-Schrodinger perturbation series; rather, they attempted to 
/it the small-e and large-a; behaviors of iplx) and then to deduce formulas for the eigenvalues. 
Their work represented the wave function iIj{x) by an exponential function only, whereas 
in our work there are always polynomials V multiplying the exponential structure. These 
polynomials will never disappear, even after an infinite number of resummations, because 
they express the physical- optics correction to the geometrical-optics approximation to the 
wave function. The exponential in Eq. (|4.3|) is only the geometrical-optics approximation. 



V. MULTIPLE-SCALE PERTURBATION THEORY APPLIED TO THE 
SCHRODINGER EQUATION FOR THE ANHARMONIC OSCILLATOR 



In this section we attempt to obtain the large- a; asymptotic behavior of the wave function 



^p{x) by applying MSPT directly to the Schrodinger equation ( |4.1| ). We will see that multiple- 
scale analysis gives the results of the previous section but it avoids the need to resum the 
perturbation series. 
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Recall that in Sec. |^ we showed that the asymptotic form of the wave function as the 
exponential of a cubic could only be obtained after an infinite number of resummations. Since 
multiple-scale analysis is equivalent to resumming the most secular terms in the perturbation 
series, we expect that at every order of MSPT we will find a result similar in structure to 
the forms of the previous section. Thus, we expect that leading-order MSPT based on just 
two scales will fail to generate the correct cubic asymptotic behavior of the wave-function 
in Eq. (^4.3| ). However, we expect that as we include more and more scales in the MSPT we 
will obtain a series like that in Eqs. ( [4.32| ) and (|4.33|) . 

We will see that the application of MSPT to the Scrodinger equation is not straight- 
forward. MSPT is conventionally applied to systems that reduce to classical harmonic oscil- 
lators when the perturbation parameter vanishes so that long scales are equal to the short 
scale multiplied by increasing powers of the perturbation parameter e. The multiple- scale 
analysis that we perform here is unusual because, as we have already see in Sec. |V[ the 
long scales are now proportional to increasing powers of x"^ as well as e. In the following 
generalization of MSPT we will show how the definition of the long scales relative to the 
short scale x can be deduced from the formalism and that secular ity has a natural analogue 
in the behavior of the wave function for large x. 

To illustrate our procedure we begin by performing a multiple- scale analysis involving 
only two scales. (A higher-order analysis will be given later on.) We assume a perturbative 
solution to Eq. (|4.1j ) of the form 



where the short scale is x and the long scale is an unknown function of x: C,i = e/i(x). 
(We use the subscript 1 to indicate that this is the first of a hierarchy of longer and longer 
scales.) This assumption is reminiscent of the method of stretched coordinates 0. 
Using the chain rule we calculate the first two derivatives of ^/'(x): 



'^Pix) = Go{x,^,) + eG,{x,^^) + Oie^), 



(5.1) 
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dx'^ dx'^ \dS,idx dx d^i dx"^ dx"^ j 
Next, we substitute Eqs. ( |0| ) and ( |5.2| ) into the Schrodinger equation and collect co- 



efficients of powers of e. To order we obtain the Schrodinger equation for the quantum 
harmonic oscillator: 

d^Go fl o 1 



+ -a;' - - Go = 0. (5.3) 



dx'^ V4 2 
The solution to this equation that is normalizable is 

Go(x,el) = A(el)e-i^^ (5.4) 

where A{^i) is an unknown function of the long scale 
To order we obtain 



dx"^ \A 2/ d^idx dx d^i dx'^ 4 

where i?i is the ffist-order correction to the energy eigenvalue 
Substituting the solution for Go into Eq. ( |5.5| ) gives 



d'Gi /I 2 1 



1 ^2 



d^i \ dx dx'^ I v4 



(5.6) 



^4 2, 

We now make the crucial argument in this procedure. By analogy with conventional 
MSPT we do not want Gi to satisfy an inhomogeneous equation driven by its homogeneous 
solution Go because this will give a contribution to Gi that may decay less quickly than 
Go for large x. For example, if the square bracket in Eq. (|5.6| ) were a constant, then Gi 
would grow like Go log x for large x. We thus require that the square bracket on the right 
side of Eq. ( |5.6| ) vanish identically. Of course, one cannot give a logical argument that 
this expression must vanish. Rather, we exploit the functional degree of freedom that we 
introduced in Eq. ( p.l| ) when we replaced a function of one variable by a function of two 
variables and simply demand that the square bracket vanish identically. 

We proceed as one does when performing a separation of variables for a partial differential 
equation. The vanishing of the expression in square brackets gives two ordinary differential 
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equations, one for /i and one for A. Apart from constants of integration that are irrelevant 
because they scale out, we obtain 

fJx) = —x^ + -x' and ^1 = -, (5.7) 

which fully defines the scale ^1 = e (^j^x"^ + fx^j. Note that Ei computed here is the usual 
perturbative correction to the eigenvalue determined from Eq. ( [4 .81 ). Also, we have 

A(ei)=e-«\ (5.8) 

Finally, we use the expression for ,^1 in terms of x to rewrite the solution for Go'- 

Go = e-i^'(i+i^)-i^^^'. (5.9) 

Observe that this generalized multiple-scale procedure leads to a solution that incor- 



porates not only the next term in the expansion of the exponent [see Eq. (|4.32| )] but also 
contains a correction to the term. If we expand the exponent, we see that this correction 
gives the constant term in the polynomial Vi^\zi^) in Eq. ( 4.45| ). 



We now perform the multiple-scale approximation to higher order. We will see that our 
multiple-scale procedure determines the new long-scale variables required to describe the 
behavior of the wave function for large x. Moreover, the higher-order calculation predicts 
exactly the higher-order corrections to the energy eigenvalue. We generalize the first-order 
calculation by introducing more long-scale variables ^2 = ^^f2{x), ^3 = e'^/3(x), and so on. 
The functions fn are determined from the requirement that the right side of the differential 
equation for G„ vanish. This requirement yields the new dependence of Go on 

To order this procedure yields 

2 / 1 6 11 4 21 - - 21 



and 



= i—x^ + —x^ + —x'] and £2 = , (5.10) 

^ \96 64 8 / ^ 8 ' ^ ' 



Go(x) = e-^^H^+i^-t^^)-l^-1l-x^)+^^^^^ (5.11) 

To order we have 
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s / 1 8 21 fj 45 4 333 2\ , ^ 333 

.3 J -r6 I 4 , 2 



£3 = eM + + — + and = , (5.12) 

V256 192 32 32 y ^ 16 ' ^ ^ 

and 

(^0(0;) =6 4-^ [,'-^2'^ 4 + 8 J le*^-^ 4''^ 2'^ J^96^ ^ 2 "^J 256 *^ ^ . (5.13) 



Observe the pattern that develops. Every new order in MSPT reproduces a new term in the 
expansion of the exponent in Eq. ([4.32|) together with x-dependent corrections to the wave 
function in the exponentiaL To recover the sequence polynomials Vi, V2, V3, . . ., which are 
discussed in Sec. |V| one need only expand the exponentials. 
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